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I. INTRODUCTION 


A. BACKGROUND 

Personnel assignment problems in the armed forces usually 
have two major objectives: Fill as many vacant positions as 
possible, and fill the positions with the most qualified 
persons. Different positions require different qualifications, 
and because of more and more sophisticated weapons and 
equipment, a large percentage of military personnel needs 
occupational training to acquire these qualifications. To 
minimize training time and cost while furthering the second 
objective, individuals must be assigned to positions in a 
manner which makes optimal use of their existing and potential 
abilities. 

In this connection the initial assignment of recruits to 
occupational training is of special importance. It is not only 
to be seen under the above mentioned objectives but also with 
respect to its influence on individual careers and job 
satisfaction by 'pushing the recruit in a certain direction’. 
Therefore, the decision process determining the initial 
assignment must not be based on rule-of-thumb procedures 
and/or intuition but on effective alternatives resulting from 


a profound problem analysis. 


The Manpower Management Enlisted Assignment (MMEA) Branch 
of Headquarters United States Marine Corps (HQMC) faces this 
decision process about once a week, when on average 750 
recruits are to be distributed among 60 or more training 


classes after graduation from the School of Infantry (SOI). 


B. GOALS OF THE THESIS 

The purpose of this thesis is to develop a prototype 
computer system based on a mathematical optimization model 
which can be used to help MMEA assign Marine recruits to 
occupational training. It must have the following properties: 

1. The model should give realistic answers which can be 
directly usable for nearly all recruits considered. 

2. The model should be as generic as possible. An increase 
in problem size (e.g. more recruits) must not result in 
any model changes or in extremely degraded performance. 

3. Only integer solutions are acceptable. 

4. The model must be easy to modify in order to accommodate 
changes for future concerns. For this reason, it should 
use general purpose solvers, rather than rely upon 


algorithms designed for a unique purpose. 


5. The results of the model must be directly importable into 
other software packages (e.g. spreadsheets). 


6. The implementation should not be hardware specific. The 
model should run on mainframes, workstations, and 
386/486 Personal Computers, and it should also be 
implementable in future computer environments. 


C. TERMINOLOGY 

This section describes terminology used in the model. 

1. Prerequisites 

When recruits are selected for training classes, their 

existing qualifications are compared with desired 
qualifications which are defined for each class. The desired 
qualifications are called prerequisites. We distinguish 
between mandatory and desirable prerequisites: Mandatory 
prerequisites describe the minimum qualifications which a 
recruit must meet to be eligible for a training class, while 
desirable prerequisites define qualifications which a recruit 
should have in addition to the mandatory prerequisites, in 
order to be a more desired candidate for a training class. 

2. Prerequisite Level 

For each training class the prerequisites are combined 

in one or more sets, where each set represents a level of 
desired qualification of the trainees. The sets are called 
prerequisite levels. The levels build a hierarchy of quality 
and are ranked by integers, starting w ,1 one for the set of 
mandatory prerequisites as the lowest level. In order to 
satisfy prerequisite level n of a class a recruit must have 
not only the qualifications described in level n, but also 
the qualifications of all lower levels n-1, n-2, ..., 1. A 


recruit who is not eligible for a class is said to have 


prerequisite level zero for this class. The number of 
specified prerequisite levels can differ between classes. 
3. Quota 
The number of seats in a training class is called a 
quota. The quota is the largest and at the same time the most 
desirable number of trainees to enroll ina class. 
4. Fill Priority 
The fill priority of a training class is an integer 
which ranks the relative importance of filling the class. Fill 
priorities are necessary inputs because of the frequently 
occurring possibility that not all quotas can be met. Classes 
with the same fill priority form a priority group. 
5. Fit Priority 
The fit priority of a training class is an integer 
which ranks the relative importance of filling the class with 
recruits having high prerequisite levels. 
6. Area Aptitude Composite 
Each Marine recruit must pass the Army Area Aptitude 
Battery in which his or her intellectual, psychic, and motor 
abilities are tested. The Area Aptitude Composites (AAC) are 
linear combinations of the scores a recruit received in the 
basic tests of the battery. The AACs are used to estimate a 
recruit's success in the training classes. For each class one 


relevant AAC is specified. 


D. THESIS OUTLINE 

In Chapter II, the Marine Corps objectives and policies 
for recruit assignments are presented. The current solution 
and data sources are also described in this chapter. 

A basic recruit assignment model is formulated in Chapter 
TII. The difficulties in finding feasible, integer solutions 
to this model are also discussed. The basic model is not 
presented as a viable approach, but as a valuable initial 
framework. Chapter IV contains the redevelopment of the basic 
model into a two phase model. The recruit assignment problem 
is decomposed into a linear programming subproblem and a 
network subproblem , which are solved in sequence. 

The computer implementation of the model, preprocessing of 
input data, and results are described in Chapter V. 
Conclusions and recommendations for future improvement are 


given in Chapter VI. 


II. ASSIGNMENT OF MARINE RECRUITS TO OCCUPATIONAL TRAINING 


A. OBJECTIVES OF HEADQUARTERS USMC 

The Manpower Management Enlisted Assignment Branch at 
Marine Headquarters has four objectives for the recruit 
assignment problem: 

1. Maximize the fill of the training classes. If quotas 
cannot be met, allocate vacancies according to the fill 
priorities. 

2. Maximize the quality of assignments as measured by 
prerequisite level, with competition for high quality 
recruits adjudicated by the fit priorities. 

3. Minimize the total waiting time between the recruits' 
graduation from Infantry School and the beginning of 
training classes. 

4. Maximize the expected success of the recruits by 
assigning each recruit as close as possible to the 
training class for which he has the highest proficiency 
as measured by the Area Aptitude Composite (AAC). 

The objectives are listed above in their current order of 
importance. The models developed in this thesis would allow 
for the order to change. 

Objectives 2 and 4 are often correlated, since for all 
training classes a minimum score in the relevant AAC is a 
mandatory prerequisite, and for most training classes the 
required score increases with the prerequisite level. On the 


other hand, these objectives can also be conflicting. A 


recruit might not be eligible for the training classes 


corresponding to his or her highest AAC, because he or she 

does not meet other mandatory prerequisites of the classes. 
The waiting time objective depends on the assignment 

decisions because the recruits can have different availability 


dates and the training classes have different start dates. 


B. POLICIES OF HEADQUARTERS USMC 

Besides the constraints that each recruit has to be 
assigned to exactly one training class and that each recruit 
must meet the mandatory prerequisites of his or her assigned 
class, three policies have to be observed in the recruit 
assignment problem. 

1. Program-Enlisted-For (PEF) 

A high percentage of Marine recruits are guaranteed by 
their recruiters that they will get a job out of a specific 
group of assignments. These agreements must be honored. The 
given PEF guarantee can possibly overrule the constraint on 
mandatory prerequisites for training classes. 

2. Minority Distribution Policy 

Recruits are classified as either minority group 
members or non-minorities. Each training class must adhere to 
a specified minimum and maximum percentage filled with 
minorities. The percentages can differ between training 


classes and can also change with time. 


3. Allocation of Shortages 

While the f1i11 priorities guide the allocation of 
shortage of qualified recruits between priority groups, 
sharing coefficients determine the allocation of shortage 
between training classes within the same priority group. For 
each class j a sharing coefficient s, > 0 specifies desired 
sharing targets within priority groups as follows: 

Let classes k and j belong to the same priority group. 
If all quotas in this priority group cannot be fully 
satisfied, a desired outcome is to have classes k and j share 
the shortage so that the ratio of their shortage percentages 
is equal to the ratio of their sharing coefficients, i.e. 
( shortage, %* ) / ( shortage, % ) = s / s;.- Typically, the 
sharing coefficients are one, so that equal sharing of 


shortage within priority classes is encouraged. 


C. CURRENT SOLUTION 

MMEA currently uses the "Recruit Distribution Model (RDM)" 
which is a product of the Decision Systems Associates, Inc. 
(DSAI) [Ref. 1]. The model was first developed in 1965 and has 
been adapted by DSAI to changing objectives and policies since 
then. RDM consists of a system of optimizing algorithms, which 
are applied sequentially to the problem, generating at each 


stage additional constraints for the following stages. RDM 


currently runs only on a Control Data Corporation CYBER 175, 
a mainframe computer that was introduced in the early 1970s. 

The underlying mathematical model of RDM has not changed. 
As new concerns of the Marine Corps were encountered, many of 
them were incorporated through the definition of prerequisite 


levels. 


D. DATA SOURCES 
1. Recruit Data 
Recruit data records are kept in the USMC Recruit 
Accession Management System (RAMS). These records contain the 
name, Social Security number and all characteristics of the 
recruit which are necessary to determine his or her 
eligibility for each training class. 
2. Training Class Data 
Quota, start date, and fill priority of each training 
class are provided in the Training Quota Memorandum (TQM), 
which is generated by Marine Corps Development and Education 
Command. Fit priorities, codes of relevant AACs, number of 
prerequisite levels, and the Sys set for each 
prerequisite level are maintained by MMEA ina file called the 
Prerequisite Dictionary. 
3.  Program-Enlisted-For Data 
MMEA maintains a catalog of all PEF guarantees and the 


corresponding promised assignments. From this file we obtain 


the information to which training classes a recruit with PEF 


guarantee can be assigned, and for which of these classes the 


constraint on mandatory prerequisites may be overruled. 
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III. BASIC RECRUIT ASSIGNMENT MODEL 


Pes FORMULATION OF THE MODEL 

This section describes an initial translation of the 
recruit assignment problem into a mathematical model. This 
model is not regarded as a viable approach on its own, but it 


serves as a framework upon which a usable model can be built. 


Indices 
i1eéetr - available recruits 
Jed - training classes 
Sets 
GC: - set of all classes for which recruit i 
is eligible 
R, - set of all recruits who are eligible for 
class j 
MI - recruits who are members of a minority 
MIR, - intersection of MI and R; 
Given Data 
q; - quota of class j} ( men ) 
lmin, - minimum minority fraction of class j 
hmin, - maximum minority fraction of class j 


I 


S, - sharing coefficient of class j 


filpr, - fill priority of class j 

fitpr, - fit priority of class j 

aac,; - score of recruit i in AAC relevant for 
class j 

ti; - time between graduation of recruit i and 


start of class j 
pl,; - prerequisite level which recruit i meets 


for class j 


Decision Variables 


Xi 1 if recruit i is assigned to class j 


O otherwise 


Constraints 


Yo xi; = 1 for all iel (3 


JE Cy 
Assign each recruit to exactly one class. 


yo x aes for all jeJ (Jaze 


Observe the upper limit of seats (quota) for each class. 


12 


sibel a < hmin, for allje gd (3.3) 


Le “3 
Le Ry 


Observe the minority distribution for each class. 


qj = ij 
i€ Ry 
S, 
25) ae ed forval bear ikea (334) 
Sx 
Qk 7 23 Xik 
Le Ky 
or 


where classes j} and k belong to the same priority 
group and k < j. 
Observe the policy for share of shortage between classes 
with the same fill priority. 


Constraints (3.3) and (3.4) can be rewritten in linear 


form: 
Imin, > x; < > xi; < hmin,; )O x; Gea) 
6 K, 1€ MIR; 1€é Ry 
qj - x Qk ~ Xik 
4 i€ R; Ee ié€ R, : (3.4') 
qd; on 
pen(Ol) pitor all (ig) 6(1,0)) te.6o 16, (3.5) 


The decision variables are binary variables, and are 


defined only if recruit i is eligible for class j. 


r3 


Objective Functions 


The problem has four objective functions: 


Maximize > filw, > x 


jeg 1€é Ry; 


Maximize fill of training classes according to the fill 


DEIOrULies. 
Maximize fitw; }> pl,; x,; (327) 
Jed PEeR, 
Maximize quality of assignments according to fit 


priorities and prerequisite levels. 


Minimize YY Deer, (3.8) 
JO J LEK, 


Minimize waiting time. 


— aac.. 
Maximize yy — «x,, (3.9) 


ij 


Maximize expected success of recruits. 


Parameter maxaac, is the highest score recruit i achieved 
in all AACS (maxaac,; = max {aac,, | Jj € J}). 

Parameter filw, is the fill weight of class j and is 
calculated by: filw, = ((minfilpr + 1) - filpr)*, where 


minfilpr is the lowest fill priority. Squaring the difference 
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accentuates the requirement that classes with high fill 
priorities are filled first. 

Parameter fitw, is the fit weight of class j and is 
calculated by fitw,; = (minfitpr + 1) - fitpr,, where minfitpr 


is the lowest fit priority. 


B. DISCUSSION OF THE BASIC MODEL 
1. Objective Functions 
The recruit assignment problem is a multiobjective 
optimization problem for which various approaches exist 
(Ref. 2]. One approach is to specify a weight for each 
objective function according to its relative importance, add 
the weighted functions, and solve the problem using a solver 
for single-objective problems. This procedure will be applied 
to the basic model. The objective functions of the basic model 
have different units and differ in the size of the parameters, 
which means they are not comparable. Therefore, the objective 
functions must be transformed to a common scale before they 
can be weighted and added. 
2. Feasibility 
The model described above 1s seldom feasible. A 
necessary but not sufficient condition for feasibility is that 
the total number of recruits is less than or equal to the sum 
of the quotas. Since the quotas are often planned before 
enough information about the recruit resources is known, this 


feasibility condition can be violated. 
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Another source of infeasibility is the PEF policy, 
which heavily reduces the number of eligible recruit-class 
combinations. This can result in the situation that some 
classes have more eligible candidates than their quotas, while 
other classes have a drastic shortfall. The situation is made 
even worse by the sharing constraint (3.4), which, if 
enforced, would prevent all classes of the priority group from 
getting a higher fill percentage than the class with the least 
fill. Consequently, in practice, either some eligible recruits 
do not get an assignment or the sharing constraint is 
violated. Either case is infeasible in the initial model, but 
this situation may not be avoidable. 

Also, the minority constraint may not be satisfiable 
because the number of qualified minority members can fall 
below the requested lower limit or exceed the upper limit. 

The basic model has to be modified in a manner which 
takes the above mentioned possibilities for infeasibility into 
consideration and guarantees usable results. The model should 
choose which constraint to violate based on policy parameters 
entered by the user. 

3. Integrality 

The model is a zero-one integer program (IP). Without 
the minority and sharing constraints it is a network problen, 
for which integer solutions are guaranteed using a linear 


program (LP) solver (Ref. 3]. These two constraints destroy 
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the network structure and cause the solution of the 
corresponding LP to fractionate. This means that, in order to 
get integer solutions, either an IP algorithm or a heuristic 
must be used. The solution time of IP problems can drastically 
increase as the number of variables increases. The typical 
recruit distribution problem has about 750 recruits and 60 
training classes with about 20% of the recruits, on average, 
eligible for each class. This problem size yields a model size 
of 9000 binary variables, which is much too large to guarantee 
acceptable run times, given the model's structure. Therefore, 
a rounding heuristic applied to the solution of the LP 
relaxation seems to be favorable for this problem. We have 
chosen an optimization-based heuristic. 

In the next chapter we develop a new, practical 
approach which overcomes the shortcomings of the initial, 


basic model. 


ay, 


IV. DEVELOPMENT OF THE MODEL 


A. ELASTIC VARIABLES AND PENALTY COST 
1. Mathematical Background 

A common procedure to overcome infeasibility is to 
elasticize (sometimes called "soften") the constraints which 
can cause infeasibility [Ref. 4, 5]. Elasticizing a constraint 
means introducing additional nonnegative variables which 
represent the under- or overachievement of the originally 
desired range of the constraint. These "elastic" variables are 
multiplied by penalty costs and added to the objective 
function. 

Elasticizing constraints yields useful information 
when the original model is infeasible. Nonzero elastic 
variables in the optimal solution indicate which constraints 
cause infeasibility in the original model and, by extension, 
which parameters would have to be changed in order to move 
toward feasibility. When it is possible to - avoid 
infeasibility, the elastic penalties enable the model to 
select which constraints to violate and by how much. It also 
allows the model to reflect the common practice of trading off 
satisfaction of one constraint for another constraint or for 


objective function improvement. 
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2. Elasticizing the Minority Constraints 

The minority constraints (3.3') are elasticized by 
subtracting the positive variable dminu,, which is the 
shortage of minority members in class j, from the left hand 
side of the inequality, and adding the positive variable 
dmino,, which gives the number of minority members exceeding 
the upper minority limit of class j, to the right hand side. 
To keep the sum of the penalty costs of this constraint 
between zero and one, the penalty costs must be divided by the 
largest value that dminu, or dmino, can take on. That value is 
maxm, = max {ming, q;, (1 - hming;)q,;}. So pmin,; is specified 
for all j € J by pmin, = 1 / (maxm, A), where A is the total 
number of planned classes. 

The new constraints are: 


Imin; )) xi; - DD Xiz - dminu, < 0 for alljeJg 
LER; i€ MIR, 


(4.4) 
Observe the lower bound on minorities for each class. 


Mga me > 5a 


a € MIR, i€ R; 


ee dmino, < 0 for alljeJ 


(4.5) 
Observe the upper bound on minorities for each class. 
3. Elasticizing the Sharing Constraints 
The linear sharing constraints (3.4') are elasticized 


by adding to the left hand side the positive variables dsu,, 
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and dso,,, which represent underachievement and 
overachievement, respectively, of the desired fill percentage 
of class k in comparison with class j. Since both violations 
are equally undesirable the same normalized penalty costs ps,, 
are specified for both variables by ps, = 1 / (smax * N), 
where N is the total number of combinations (j,k) € (JxJ), so 
that classes j and k belong to the same priority group and 

k < j, and smax = max {s; | J € J}. The new constraint is: 
qd; ~ De Xij dy ~ > Xik 

i€é R; i€ Ry 


-S. ——*—- + dsu., -dso,, = 0 
q; y, ay jk jk 


(4.6) 


for all j, k € J, where class j and class k are in the 
Same priority group and k < j. 
Observe the policy for share of shortage between classes 
with the same fill priority. 
4. Elasticizing the Supply Constraints 
The supply constraints (3.1) are elasticized by adding 
to the left hand side of the equation the positive variable 
da,, which allows the possibility that recruit i is not 
assigned to any class. When determining the penalty cost pa, 
we must take into consideration that only eligible recruits 
can be assigned to a class, and that recruits with a PEF 
guarantee must get an assignment first. The following costs 
satisfy these conditions: pa, = 0, if C, is empty; pa, = 3, if 


C, is not empty and recruit i has no PEF guarantee; 
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pa, = 3 * NOPEF, if C, is not empty and recruit i has a PEF 
guarantee, where NOPEF is the number of all recruits without 
PEF guarantee. The penalty cost for not assigning eligible 
recruits is higher than the sum of all other penalty costs. 
This causes the model to assign eligible recruits, even if the 
Original minority and/or sharing constraints have to be 
violated. 
The new constraint is: 


> Xi; + da, =1 for all iel (4.2) 
sre 


Assign each recruit to at most one class. 


B. NORMALIZED AND WEIGHTED OBJECTIVE FUNCTIONS 
1. General Idea 

A widely used procedure to make conflicting 
objectives comparable is to score each of them on the scale 
(O, 1]. Let A, be the achievement level of objective n. Define 
a, = (A, -— Aw) / (As - Aw), where A, is the most desirable and 
A.w the least desirable achievement lev: ” of objective n, then 
a, € [0, 1]. This method, called "proportional scoring", will 
be applied to the objective functions of the basic model. The 
scored functions will be weighted in accordance with their 


relative importance. 


Ze 


2. Fill of Training Classes 

It 1s most desirable that the number of assignments be 
equal to the sum of the quotas, and least desirable, that no 
recruit be assigned to any class. Therefore the fill objective 
is normalized by dividing by the sum of the quotas. 

The fill weights filw, as described in Section III A 
not only determine the order in which the classes are filled, 
but also weight the fill objective function relative to the 
other objectives. 

3. Fit of Assignments 

The most desirable level of the fit function is 
achieved if all recruits in a training class have the highest 
prerequisite level of this class. The least desirable 
achievement occurs if all recruits have only the mandatory 
attributes, i.e., prerequisite level one. (Level zero is not 
possible, since those recruits are not eligible for this 
Class.) Therefore, pl,; 1s replaced by 
( pl,,; -1 ) %& ((hl, - 1)R) in function (3.7), where hl, degen 
highest prerequisite level of class j and R is the total 
number of eligible recruits. If hl, = 1, class j is omitted 
from the fit objective, because the fit quality is not 
controllable for such a class. 

The fit weights fitw, as described in Section III A 
serve two purposes. They give the classes priorities for 


receiving trainees with high prerequisite levels, and they 
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give weight to the fit objective function relative to the 
other objectives. 
4. Waiting Time 
The optimal waiting time for all assigned recruits is 
zero, while the worst case 1s maxt = max {t,,| i ¢€ I, j € J}. 
Therefore, t,; 1s replaced by -(maxt - t,;) / (maxt R) in 
munction (3.8). 
5. Expected Success of Recruits 
A recruit is most likely to succeed in the training 
class for which he or she has the highest proficiency, as 
measured by the Area Aptitude Composite. The expected success 
decreases as the AAC score decreases and is defined to be zero 


if the score is zero. For scaling, aac,, / maxaac, in function 


(3.9) 1s replaced by aac,, / (maxaac, R). 


=. FORMULATION OF THE ELASTICIZED MODEL 
The model with elasticized constraints can be written as 


follows: 
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Maximize 


’ ake 
JET i€ Ry 
2, Fj 
JET 
pl, iL 
+ Ww fitw. > J x 
fit a, _ 5 
JET ié€ R, (hl; 2) Rk 
maxt - ae 
+ W.- 
time ij 
- cae ie P rae axaac, R14 
je Ji€é Ry 


~ — Wsupp y pa; da; 
Lean 


jéed 


pe Wohar DS 5, (dsu,;, a dSO;, ) 


JE TkET 


(4.1) 


(The seven parts of the objective function are weighted as 


described in Section IV B. The parameters wy}, Writ, etc. were 


all set to one in our runs of the model, but they are included 


to give the user an additional possibility to tune the model.) 
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Subject to 


Y xX, + da, =1 for a lla es (4.2) 
FEC; 


Assign each recruit to at most one class. 


Der x (= a: for ail je (4.3) 
fear 


Observe the quota of each class. 


Torney is eg) mn, <0 for alljeJ 
Te R; Pe MIR; 


(4.4) 


Observe the lower bound on minorities for each class. 


x,,; - hmin, > x,; - dmino, < 0 for all jéa 
i € MIR, GER, 


(4.5) 


Observe the upper bound on minorities for each class. 


2S Xij fe es 
(4.6) 
Ss cas, Sg, — Ee a dai. =ds0, = 0 
k = j a Ae Po 
| 


for all j, k € J, where class j and class k are in the 
same priority group and k < j. 
Observe the policy for share of shortage between classes 


with the same fill priority. 


da. 20 for alliel 


22 


dminu,, dmino,, dsu;,, dso;, 2 0 for allj,ked 


xy © (Oye for all i1é1, j) € J,)s.t. ig 


D. ‘'TWO-PHASE APPROACH 
1. Motivation 

The elasticized model is still an IP problem for which 
an LP solver is very likely to yield fractional solutions. To 
get completely integer solutions, the following two-phase 
approach is taken: 

In the first phase, the LP relaxation of the 
elasticized model is solved. This means the integrality 
constraint is replaced by x,;; 1s greater than or equal to zero. 
The optimal solution of the LP relaxation is used to calculate 
integer upper and lower bounds’ on the fill of the classes. 
These bounds are kept as close as possible to the relaxed 
problem's optimal fill values, so that, as long as the fills 
are varied only within the bounds, the corresponding minority 
and sharing ratios will also stay very close to the values 
found by the optimal solution of the LP relaxation. 

Therefore, in phase two we can replace the sharing 
constraint by using the new bounds in the fill and minority 
constraint. This results in a network model, which guarantees 
completely integer LP solutions. 

A similar two-phase approach was used quiet 


successfully in scheduling flowlines [{Ref. 6]. 
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2. Computation of the Bounds 


Let 


fill, = ee erect For alles teed 


where x,,‘"? are the optimal solutions of the LP version of the 
elasticized model. The lower bound of the fill of class j is 
then defined by log; = floor (f1i11,), where floor (a) is the 
greatest integer smaller than or equal to a. The upper bound 
is defined by hiq, = min {loq, + 1, q,;}, which makes sure that 
the final fill value is no more than one unit away from phase 
one fill, and also no greater than the original quota. 

Based on these bounds the smallest and largest number 
of minority recruits is calculated. The lower limit is given 
by lom, = floor (lmin, hiq,;), and the upper limit by 
him, = ceiling (hmin, hiq,), where ceiling (a) is defined as 
the smallest integer greater than or equal to a. 

3. The Network Model 

Mathematically the network mc~7l can be stated as 

follows: 
Maximize 


»s yy Uig *i; > Weupp Dee da ; (4.7) 
Le 


JETIER; 


- Wain > pmin,; (dminu,; + dmino, ) 
fev 


Pad | 


where u,, 1s the sum of the coefficients of the first 


four parts of the objective function (4.1). 


Subject to 


Yo x, + da, =1 for alliel (4.8) 
JE Cy 


Assign each recruit to at most one class. 


ee) = ~6Log: for allj eg (4.9) 
i€é R, 


Observe the lower bound on quota for each class. 


Y> xi, < hig, foralljedg (4.10) 


i€é R; 


Observe the upper bound on quota for each class. 


X;; + dminu; 2 om, for allj € J (4.11) 
i € MIR; 


Observe the lower bound on minorities for each class. 


Xi, — GMinon  < Vian for all j ES 3G) 
i € MIR, 


Observe the upper bound on minorities for each class. 
All variables are nonnegative. 
Figure 1 illustrates a network model which is 
equivalent to the phase two linear program. The nodes P, 


represent the available recruits, M,; the minority quotas, and 
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Q, the quotas of the training classes. The final node § is 
used in the network solution process to draw recruits through 
the system. Without loss of generality we can assume that 
nodes P, through P, represent minority members while P43) 
through P, denote non-minorities. Each recruit node has a 
supply of one, indicating that each recruit can be assigned to 
only one class. Node S has demand R which is the number of all 
available recruits. All other nodes have supply equal zero and 
serve as transshipment nodes. 

An arc representing x,, exists if recruit i is eligible 
for training class j}. This arc has tail P,, and if recruit i 
is a minority member, it has head M,. If the recruit is a 
nonminority, it has head Q;. This arc has objective function 
coefficient u,,, the gain of assigning recruit 1 to class j. 

The arc from M, to Q, represents the required minority 
enrollment in class j. The lower and upper bounds on this arc 
are obtained from the solution of the phase one subproblem as 
described above. Similarly, the arc from Q,; to S represents 
the required total enrollment in class j. The arcs (M,,Q;) and 
(Q,;,S) have zero cost because traversing them means that 
constraints are being satisfied. 

The elastic variable dmino,, corresponding to excess 
minorities in class j, is represented by the arcs from 
minority recruits to Q,. These arcs have an objective 


coefficient of -pmin, representing the elastic penalty cost 


Zo 
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Network Model for Phase Two 


Figure 1. 
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for excess minorities. The possibility of a shortage of 
minorities 1s modeled with arcs from the nonminority recruit 
nodes to M,. If this elastic are is used, a nonminority 
recruit will be counted in a minority quota, for which a 
penalty of pmin, is charged. 

Each recruit is also connected by an arc to node S, so 
that if recruit i cannot be assigned to any class, he will be 
assigned directly to node S at penalty cost pa,. 

4. Summary of the Two-Phase Model 

The two-phase model is solved as follows: 

1. Solve the LP relaxation of the elasticized model defined 
by (421) through (4.6). 


2. Compute the lower and upper bounds on minority fill and 
total fill by the formulas in Subsection IV C 2. 


3. Solve the network model defined by (4.7) through (4.12). 


This combination of steps is guaranteed to achieve an 


integer solution. 


ok 


A. 


V. COMPUTER IMPLEMENTATION AND RESULTS 


GAMS MODELING LANGUAGE 


A GAMS [Ref. 7) computer program was written to obtain 


optimal solutions for the two phase recruit assignment 


problem. The GAMS modeling language was used for several 


reasons: 


lie. 


Since the GAMS model representation is the same as the 
mathematical representation, the computer program is easy 
to read and understand, which is especially important if 
the users change frequently as is common in the military. 


An optimization model implemented in GAMS is very easy to 
modify. This is important in the case the Marine Corps 
decides to add new constraints in the future. For 
example, gender distribution limits can be handled the 
same way as the minority constraints. For another 
example, a desire to spread quality recruits among 
classes can be handled the same way as the sharing 
constraints. 


GAMS is an extremely convenient way to execute a sequence 
of optimization models in which the inputs to the second 
model depend on the optimal solution of the first model, 
as required by our two-phase method. 


GAMS allows the model to be formulated and maintained 
independently from the data it uses. This means that data 
can be changed or the problem can increase in size 
without resulting in problem changes. 


Due to the system design of GAMS, future solution 
algorithms can be used without the requirement to change 
the representation of the model, e.g., a faster linear 
program or network solver can be substituted for the 
current solver by changing one line of GAMS input. 


The GAMS system does not require special input or output 


procedures. Input files can be written with any text 
editor or can be the result of other software packages. 
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In the recent version of GAMS, version 2.25, the GAMS 
program can be linked to spreadsheets. This can be used 
in our problem to write the recruits' assignment orders 
using the optimal assignments found by GAMS. 


GAMS programs are portable from mainframes to 
workstations and PCs and vise versa without any changes 
in the program. GAMS is not designed for a special 
hardware. 


The above described features of GAMS satisfy the goals of 


this thesis as specified in Chapter I. 


B. 


PREPROCESSING OF INPUT DATA 


The raw data given by the data sources described in 


Chapter II are preprocessed as follows: 


IE 


Pie 


From the recruit data the scores of the Area Aptitude 
Composites, minority membership, graduation date, anda 
PEF flag are written to the file RECDATA. 


The qualifications of the recruits as given in the 
recruit data are compared with the prerequisites of the 
available training classes, which are contained in the 
Prerequisite Dictionary. By that we find the prerequisite 
level of each recruit for each class. For recruits with 
PEF guarantees the prerequisite levels of inappropriate 
classes are set to zero. The prerequisite level matrix 
(rows = recruits, columns = classes) is written to file 
PRELEV. 


Quota, start date, and fill pric .ty of each training 
class, as given in the Training Quota Memorandum, are 
combined with the code of the relevant AAC, fit priority, 
number of specified prerequisite levels, and the lower 
and upper limit of the minority fraction in file QUOTAS. 


The first two input files are generated by a FORTRAN77 


program which was written by the author, whereas the QUOTAS 


file was written with a text editor. A generic program which 


33 


generates the three input files would have exceeded the scope 
of this thesis. Figure 2 shows which raw data (rectangles) are 


used to generate the input data (ellipses) for the program. 


Catalog of 


Prerequisite 


Recruit Data 
T QM 
PEF 





Dictionary 


x 


RAMS 





—_ Ce | (om 


Figure 2. Preprocessing of Input Data 


Ideally, the raw data should reside in database management 


system (DBMS) and the GAMS program should extract appropriate 


Gata and execute via DBMS commands. Such an environment is 


described in "TC-EXPRESS Version of the Recruiter Allocation 


Model" [Ref. 8]. 
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C. COMPUTATIONAL RESULTS 
The model was tested with real data of 461 recruits and 65 
training classes obtained from MMEA. The training classes had 
23 different fill priorities and 9 fit priorities. The 
resulting GAMS model had 886 equations and 5858 variables in 
phase one (elasticized model), and 848 equations and 5587 
variables in phase two. The relatively small number of 
variables, compared to 461 * 65 = 29965 theoretically possible 
decision variables, is due to the PEF policy which reduced the 
number of eligible recruit-class pairings to 4869 (density = 
16.3%). This can be fully exploited in the GAMS model. 
The GAMS program was run on three computers: 
1. AMDAHL 5990-500 mainframe using the CP/CMS timesharing 
system with GAMS 2.19 and ZOOM 2.1. 
2. 486-based PC under DOS with GAMS 2.25 and MINOSS5. 
3. IBM RS6000 Model 530 under UNIX with GAMS 2.25 and 
MINOSS5. 
About four megabytes RAM were needed. Run times are shown 


an Table 1. 


D. COMPARISON WITH CURRENT SYSTEM 

For the data set described above the optimal solution of 
the Recruit Distribution Model (RDM) was obtained from MMEA. 
Since no complete model formulation of RDM was available to 
the author, only a limited comparison of our results with 


those of RDM is possible. We compared the achievements of the 
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four objectives specified by MMEA, with the results given in 


Table 2. 


TABLE 1. RUN TIMES OF THE GAMS PROGRAM 

















Computer AMDAHL 5990 | IBM RS6000 486 PC 





GAMS Version GAMS 2.19 GAMS 2.25 GAMS 2.25 






Solver ZOOM 2.1 MINOS5 MINOS5 ) 


Generation/Report 47 sec 61 sec 144 sec 


TABLE 2. COMPARISON OF OBJECTIVE ACHIEVEMENTS 


Jassignnente [a 


| Waiting Time 12,376 man days 12,385 man days 


We can conclude that there are no significant differences 






between the solutions of the two models. 
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In both solutions 460 recruits were assigned to training 
classes. (One of the 461 given recruits was not eligible for 
any training class.) 328 recruits (71.3%) were assigned to the 
same training classes in both solutions. For about 60% of the 
remaining 132 recruits the assignments to classes were changed 
pairwise (e.g., in our solution recruit i is assigned to class 
j and recruit 1 to class k, while RDM assigns recruit i to 
class k and recruit 1 to class j) without any change in the 
objective function value. 

Out of 23 priority groups 15 were filled with the same 
number of recruits in both solutions, six groups had a fill 
adifference of one, and two groups had a fill difference of 
two. No tendency was observed that one of the models achieved 
a significant better fill of groups with high priorities (see 


Appendix). 


Se) 


VI. CONCLUSIONS AND RECOMMENDATIONS 


This thesis describes a two phase optimization model for 
the assignment of Marine recruits to occupational training. In 
phase one an elasticized linear program with relaxed 
integrality constraints is solved to calculate the optimal 
fill of the training classes. These optimal values are used to 
compute integer lower and upper bounds on the fill of the 
classes for a network model which is solved in phase two, 


yielding an integer solution. 


A. CONCLUSIONS 

In Chapter V it was shown that the solutions of the two 
phase model are very similar to those of the Recruit 
Distribution Model which satisfy the requirements of the 
Marine Corps. The results look promising for using our model 
as a basis for an optimization system to assign recruits to 
training classes. The model provides not only an effective 
decision aid in the assignment process, but can also give 
feedback on the effect of the policies (constraints) by the 
elastic variables. This enables the user to modify constraint 
parameters in order to achieve more preferred results. The 
short run times of the GAMS program favor multiple runs. The 
implementation in GAMS makes the model independent of special 


hardware and allows the use of future computer environments. 
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It allows the model to be easily modified if new constraints 


are added in the future. 


B. RECOMMENDATIONS 

Before the model can be fully implemented as_ an 
operational tool, the preprocessing of the input data should 
be fully automated. In the data sources which were available 
to the author the qualifications of the recruits were 
differently coded in the recruit data records and in the 
Prerequisite Dictionary, e.g., "graduation from high school" 
is coded by "5" in the recruit data record and by "HSGRAD" in 
the Prerequisite Dictionary. This complicates the comparison 
of the data in a computer program. The coding in all data 
sources should be standardized. The raw data should reside in 
a database management system and the GAMS program should 
extract appropriate data and execute via DBMS commands. 

The output of the GAMS program should be linked to a 
software package which uses the model results to write 
assignment orders, calculate statistics, etc.. 

Important objects of further resea. h are the weights of 
the objective functions. The idea of combining objectives 
through a weighting function is equivalent to assuming a 
linear utility function [Ref. 2}. A nonlinear utility function 
could be found which more accurately reflects MMEA's 
preferences for making trade-offs between the various 


constraints and objective functions. In phase one a nonlinear 
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program with this nonlinear utility function as objective 
function can be solved as done by Harrison and Rosenthal 

(Ref. 9], while in phase two the network can still be used to 
get integer solutions. Further improvement of the results is 


possible with this enhancement. 
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APPENDIX. COMPARISON OF CLASS FILL 


TABLE 3. COMPARISON OF FILL BY PRIORITY GROUPS 


PELOCrcy FILLNPS FILLRDM Difference 
3 3 3 
4 1 1 
6 1 JE 
8 7 7 
9 5 5 
10 20 30 
11 4 4 
12 184 LS5 = 
13 30 29 1 
14 7 9 = 
16 10 TO 
17 22 22 
18 6 6 
19 5 6 =a: 
20 18 18 
AAT 19 ES 
23 18 18 
24 1 1 
2? 16 14 2 
26 7 8 =i 
z/ 30 30 
28 3 3 
29 32 32 


FILLNPS = fill of priority groups in optimal solution of two 


phase model. 


FILLRDM fill of priority groups in optimal solution of RDM. 
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TABLE 4. COMPARISON OF FILL BY TRAINING CLASSES 


Class No. Priority Quota FILLNPS FILLRDM Difference 


28 3 1 1 1 

20 3 2 2 2 

oo 4 1 ik 1 
47 6 1 iL 1 

40 8 3 3 3 

4l 8 4 4 4 

11 9 1 1 1 

Sy 9 4 4 4 

38 10 2 2 2 

Bo 10 2 2 2 

9 10 8 8 8 

10 10 18 18 18 

6 i Bak 4 4 4, 

12 12 aS 13 13 

55 12 86 80 80 

58 1 99 Gal o2 =?) 
57 13 3 3 2 1 
a3 13 3 3 3 

60 13 3 3 3 

Sp) 13 ik 10 10 

56 a 12 11 a 

5 14 19 7 9 =2 
8 16 10 10 10 

17 17 1 1 1 

14 17 21 21 21 

L5 18 1 1 1 

16 18 5 5 5 

23 19 9 1 2 = 
21 19 10 4 4 

24 20 A 4 4 

25 20 4 4 4 

27 20 4 t 4 

26 20 6 6 6 

65 2a | 6 6 

64 21 7 6 Hi =k 
63 21 7 a 6 1 
51 23 4 3 3 

52 23 6 5 4 1 
49 23 7 5 5 

50 23 7 5 6 = 
53 24 2 i: 1 
42 25 1 1 1 

44 25 2 2 2 
43 PES, 2 2 2 

54 25 4 
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Class No. Priority Quota FILLNPS FILLRDM Difference 


45 25 5 5 Le. 

46 25 5 5 5 

48 25 8 uf 1 

36 z© 2 x = 
309 26 2 2 2 

22 26 5 e. 5 

18 Ze A a 1 

33 27 1 ik iu 

34 27 ay a 1 

Ue 27 2 2 2 

30 Pa | 2 2 2 

32 27 3 3 3 

62 Zz. 4 4 3 uu 
cal 27 5 5 5 

61 2/7 14 a 2 —1 
a 28 6 5 3 

3 29 2 2 2 

2 29 4 4 4 

dE PAS 26 26 26 

4 30 3 

FILLNPS = fill of training class in optimal solution of two 


phase model. 


FILLRDM 


fill of training class in optimal sulution of RDM. 
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